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ABSTRACT 

Minor merging has been postulated as the most likely evolutionary path to produce the increase 
in size and mass observed in the massive galaxies since z^2. In this Letter, we test directly this 
hypothesis comparing the population of satellites around massive galaxies in cosmological simulations 
versus the observations. We use state-of-the-art, publically available, Millennium I and II simulations 
and the associated semi-analytical galaxy catalogues to explore the time evolution of the fraction of 
massive galaxies that have satellites, the number of satellites per galaxy, the projected distance at 
which the satellite locate from the host galaxy, and the mass ratio between the host galaxies and their 
satellites. The three virtual galaxy catalogues considered here, overproduce the fraction of galaxies 
with satellites by a factor ranging between 1.5 and 6 depending on the epoch, whereas the mean 
projected distance and ratio of the satellite mass over host mass are in closer agreement with data. 
The larger pull of satellites in the semi-analytical samples could suggest that the size evolution found 
in previous hydrodynamical simulations is an artifact due to the larger number of infalling satellites 
compared to the real Universe. These results advise to revise the physical ingredients implemented in 
the semi-analytical models in order to reconcile the observed and computed fraction of galaxies with 
satellites, and eventually, it would leave some room to other mechanisms explaining the galaxy size 
growth not related to the minor merging. 

Subject headings: dark matter — galaxies: halos — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 

^sO ' Accretion of minor satellites has been postulated as 
\ the most likely mechanism to explain the significant size 
! evolution (Daddi et al. 2005; Trujillo et al. 2006) of 
C") 1 the massive galaxies over cosmic time. This idea fits 
\ with many observational indirect evidence: the progres- 
, sive growth of the wings of the profiles of the massive 
fNj ■ galaxies with time (Bezanson et al. 2009; Hopkins et al. 
7-H ; 2009, van Dokkum et al. 2010, Carrasco et al. 2010), the 
(~1.5) larger velocity dispersion of the massive galaxies 
. £h ' at high-z compared to present-day equally massive ob- 
^ ' jects (e.g. Cenarro & Trujillo 2009, Cappellari et al. 

2009, Onodera et al. 2010, Newman et al. 2010, van de 
C3 • Sande et al. 2011), and the expected mass growth by a 
factor of two of the massive galaxies with time suggested 
by the observations (see e.g. Trujillo et al. 2011). 

On the theoretical side, Naab et al. (2009) conducted a 
pioneer work on exploring the minor merging effect on a 
cosmologically motivated evolution of a massive galaxy 
since z~3. Their study supported the idea that minor 
merging can explain simultaneously the size and velocity 
dispersion evolution, as well as a moderate increase of 
the stellar mass by a factor of two. This work has been 
later confirmed by many other authors (e.g. Sommer- 
Larsen & Toft 2010, Feldmann et al. 2010, Oser et al. 
2012). Nonetheless, although both the observational and 
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theoretical side seem to converge in a unique solution to 
the problem of the size evolution of the massive galaxies, 
a direct confrontation of both the theory and the obser- 
vations has not yet been conducted. Inferring the merger 
rate (e.g. Lopez-Sanjuan et al. 2011) from the observa- 
tions is not straightforward due to the large uncertainties 
in the determination of the merging time scales. Alter- 
natively, an immediate way of testing the model with 
the data is counting the number of satellites around the 
massive galaxies. These objects will likely be involved 
in the growth of the massive objects. Observations have 
matured enough to allow a robust determination of the 
fraction of massive galaxies with satellites at different 
redshifts and for different mass ratios (e.g. Newman et 
al. 2012, Marmol-Queralto et al. 2012, Lopez-San Juan 
et al. 2012). On the other hand, present cosmological 
simulations are large enough to permit the estimation of 
this fraction with accuracy. 

In this Letter, we explore the model predictions about 
the properties of the satellites surrounding massive galax- 
ies and their evolution with cosmic time. In particular, 
we compare the changes with redshift of the fraction of 
satellites, radial projected distances and stellar mass ra- 
tio of three different semi-analytical models (Bower et al. 
2006, DeLucia et al. 2007, Guo et al. 2011) based on the 
Millennium simulation with the data taken by Marmol- 
Queralto et al. (2012). We find that although the general 
trends with redshift are reproduced (i.e. the fraction of 
massive galaxies with satellites only changes moderately 
with cosmic time in the simulations) , all the models tend 



2 



V. Quilis & I. Trujillo 



to overpredict this fraction. 

2. GALAXY CATALOGUES 

We use the public releas e of two very large N -body 
simulations Mi llennium I (|Springel et al. I 120051 ) and 
Millennium-II (jBovlan-Kolchin et al. 1 120091 ). The cos- 
mology assumed in both simulations is a ACDM with 
the following parameters: O m = 0.25, Slf, = 0.045, 
n A = 0.75, n = 1, cr 8 = 0.9, H Q = 73 ferns -1 Mpc -1 . 
The two simulations use the same number of particles, 
2160 3 , but cover different volumes and, therefore, they 
have different numerical resolutions. Thus, the compu- 
tational boxes have sides of 685 Mpc and 137 Mpc, and 
particles masses of 1.18 x 10 9 M Q and 9.45 x 10 6 M Q for 
the Millennium I and II simulations, respectively. 

The dark matter merger trees are extracted from 
the simulations u sing a com bi natio n of friends- 
of-frien ds (FoF) (IDavis et al. I fl985h and SUB- 
FIND (|Springel et al. 1 1200 ID halo finders. The dark 
matter halos are converted into galaxies according to 
different semi-analytical models which can differ in 
the phenomenological recipes to introduce the gas and 
stellar components in such haloes. In this letter, we 
consider three semi-analyt ical models available in the 
Millennium database web ()Lemson et al. If2006h . 

The fi rst one is a version o f the Durham semi-analytical 
model (|Bower et al. I I2006D that implements feedback 
due to active galactic nuclei (AGN) as a manner to 
stop or delay the formation of cooling flows. Th e other 
two models are the ones by iDe Lucia fc Blaizot I ([2007) 
and IGuo et al. I (|201lT) . Both are very similar being 
the second one an improvement of the earlier version 
that modifies or extends some of the semi-analytical 
rec ipes. The new fea tures implemented in the model 
by IGuo et aT~l (|2011| ) are: the separate evolution of 
sizes and orientations of gaseous and stellar discs, the 
size evolution of spheroids, tidal and ram-pressure strip- 
ping of satellite galaxies, and the disruption of galax- 
ies to produce intracluster light. The effect of AGN 
feedback was already inclu ded in the earlier version by 
IDe Lucia fc Blaizot ( (|2007[ ). The initial mass function 
(IMF) used to estimate the stellar masses of the sim- 
ulated gala xies were: Chabrier (2003) for IGuo et al. I 
(120111) an d IDe Lucia fc Blaizotl ([20071 ) and Kennicutt 
(1983) bv lBower et al. I (I2006D. 

The models b y iBower et al. I (|2006| ) and 
IDe Lucia fc Blaizot I (|2007| ) use the Mi llennium I 
simula tion, whereas the third model by IGuo et aT~l 
(|201lD is based on the higher resolution Millennium-II 
simulation. 

3. RESULTS 

We have generated three galaxies catalogues using the 
Millennium web application. Each catalogue corresponds 
to one different semi-analytical model as previously dis- 
cussed. For each of the catalogues, we select all the mas- 
sive galaxies - that we call host galaxies - as those objects 
with stellar masses, Mh, between 10 11 M and 10 13 M . 
For each ho st galaxy, and for the sake of com parison with 
the data by lMarmol-Queralto et all (|2012t ). we consider 
a spherical region of radius of 100 kpc (physical units). 
Within this region, we compute the number of satellites 
with stellar masses 0.1 < M s /Mh < 1 and projected dis- 
tances smaller than 100 kpc. Once the satellites fulfilling 



such mass ratio and separation conditions are identified, 
we calculate the projected distance of each satellite and 
the mass ratio respect to their host galaxy. In order 
to study different mass scales, this process is repeated 
for smaller satellites with stellar masses in the range of 
0.01 < M s /M h < 1. Possible bias effects due to the 
use of different host masses or galaxy type-dependent 
number of satellites are minimized by using the same 
mass criteria for simulated and observed galaxies and by 
comparing both samples globally, without distinguishing 
between galaxy types. 

In the Figs. 1 and 2, we present the results of the 
previous analysis on the three considered galaxy cat- 
alogues. The figures are organized in three columns, 
each one showing the results of one particular cata- 
logue. For all the columns in the Fig.l , and from top 
to bottom, we display the fraction (/) of massive galax- 
ies that have at least one satellite and the mean num- 
ber of satellites (N s ) for those host galaxies that have 
them. The average projected distance of the satellites 
(D s ) and the average ratio of the satellite stellar mass 
over the host galaxy stellar mass (M s /Mh) are presented 
in the Fig. 2 organized as in Fig.l. The full circles (tri- 
angle) stand for the simulated satellites with mass ratios 
of 0.1 < M s /M h < 1 (0.01 < M s /M h < 1). The er- 
ror bars are computed as one standard deviation The 
observational data from iMarmol-Queralto et al. I ([2012D 
are overploted as red (blue) open circles (diamonds) for 
mass ratios of 0.1 < M s /M h < 1 (0.01 < M s /M h < 1) in 
both figures. The IMF used in this work was Chabrier 
(2003). These data are representative of the results found 
for different authors (see e.g. Fig. 6 of Lopez-Sanjuan et 
al. 2012), and for this reason we use it here as a basis of 
comparison with the simulations. 

The Fig. 1 clearly shows that all the semi-analytical 
models explored in this paper systematically overpredict 
at every redshift the fraction of massive galaxies hav- 
ing satellites in their vicinity. The model expectation 
is larger than what it is currently measured by a fac- 
tor ranging between 1.5 and 6 depending on the epoch. 
This factor is independent on the mass ratio considered 
to probe the satellite population (either 1:10 or 1:100). 

The average projected distance of the satellites to the 
host galaxy in Fig. 2 is well reproduced (within the errors) 
for the case of 1:10 ratio: ~ 40 kpc. However, in the case 
of 1:100 ratio, the observed satellites are closer to the 
massive objects than what simulations predicts. Finally, 
in relation to the typical average mass of the satellites 
within each mass ratio displayed in Fig. 2, the models 
are closer to the observations, particularly in the case 
of Bower et al. (2006). Another important aspect of 
the observations which agree broadly with the models is 
the fact that all the properties of the satellites explored 
here: the fraction of massive galaxies with satellites, their 
average distance to the host object and the average mass 
ratio of the satellite depend weakly on the redshift. 

Summarizing, model predictions agree qualitatively 
well with the observations except for the fact that they 
overpredict on average the fraction of massive galaxies 
with satellites around them by a factor of two. 

4. DISCUSSION 

The most striking of our results is the overabundace 
of simulated satellites around massive galaxies by a fac- 
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Fig. 1. — The columns stand for the results of three galaxy catalogues based on different semi-analytical models. For each model, and 
from top to bottom: the fraction of massive galaxies that have at least one satellite within a sphere of lOOfcpc radius and a projected 
distance smaller than 100 kpc, and the average number of satellites per massive galaxy when they have one of such objects around. The 
full circles (triangle) stand for the satellites with stel lar mass ratios of 0.1 < M s /M;, < 1 (0.01 < Ms/M^ < 1). The error bars represent 
one standard deviation. The observational data from Marmol-Qucralto et al. (2012) are overploted as red (blue) open circles (diamonds) 
for mass ratios of 0.1 < M 3 /M h < 1.0 (0.01 < M s /M h < 1). The local observational reference (z=0.1) from Liu et al. (2011) for the 
fraction of massive galaxies with satellites with mass ratios of 0.1 < Ms/M^ < 1 is plotted as a green open triangle, no data are available 
for smaller satellites. 
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Fig. 2. — From top to bottom, the mean projected distance of the satellites to their host galaxy and the mean ratio of the satellite stellar 
mass over the host stellar mass for the three semi-analytical samples. The symbols and colors are defined as in the Fig.l. 
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tor between 1.5 and 6 compared to the observations. At 
z=0 this is quite surprising as the stellar mass function 
of the Millennium galaxies (Guo et al. 2011) is perfectly 
compatible with the observations (Perez-Gonzalez et al. 
2008). Guo et al. (2011) uses the same observational 
database than in Marmol-Queralto et al. (2012), so the 
origin of this discrepancy can not be explained as a con- 
sequence of a different dataset. It is also irrelevant the 
slightly high power spectrum normalization in Millen- 
nium (us = 0.9), as no important effects are expected at 
the scales studied in this work (see for instance Zentner & 
Bullock 2003). The only way to reconcile both results is 
considering that the spatial distribution of less massive 
galaxies is not the same in the real and in the virtual 
Universe. In particular, the satellite galaxies are more 
clustered (a factor of two on average) around the massive 
galaxies in the semi- analytical samples than what the ob- 
servations suggest. It is worthy to note that at z > the 
semi-analytical models considered here do not match the 
observed mass functions, and indeed, they overproduce 
the number of low mass galaxies. 

The fact that semi-analytical satellites appears to be 
on average a factor of two more common around the mas- 
sive galaxies than in the observations could point out to 
a factor of two larger merger time scale in the models 
than in the real Universe (a scenario already discussed 
in the literature and which have eluded so far a conclu- 
sive answer, see e.g. the discussion in Lopez-Sanjuan et 
al. 2012). Other possible sources of the discrepancy be- 
tween the abundance of real and virtual satellite galaxies 
could be related with numerical resolution effects or the 
use of an unrealistic modeling of the galaxy formation 
processes - for moderate and small masses - in the semi- 
analytical models considered in this Letter. In order to 
discard possible resolution effects, we have performed an 
extra analysis using two galaxy catalogues produced by 
the same semi-analytical model described in Guo el al. 
(2011) on both Millennium I and II simulations. The 
results are extremely similar for all the analysed quanti- 
ties in both cases. This fact discards uncontrolled effects 
due to the numerical resolution and stresses the crucial 
role of the considered semi-analytical model. In partic- 
ular, the different treatment for identifying the central 
galaxies and for the gas stripping between the models by 
Guo et al. (2011) and DeLucia et al. (2007) should en- 
hanced the satellite galaxy population in Guo's model as 
it is proved by the results shown in Figure 1. We note, 
however, that even with such differences the three mod- 
els considered here show a clear excess in the fraction of 
satellites around massive galaxies. 

Another remarkable result of the observations is the 
constancy of the fraction of massive galaxies with satel- 
lites at all redshift. This suggests that the number of 
satellites per host is in equilibrium leading to a constant 
accretion of stellar mass by the host galaxy (Nierenberg 
et al. 2012). This is in contrast with the decline of this 
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fraction since z ~ 1.2 in the higher resolution Millennium 
II simulation. We remark, however, that in the lower 
resolution Millennium I simulations this fraction seems 
also not to change with cosmic time in better agreement 
with the observations. We can link this finding with the 
radial distribution of the satellites. In all the three semi- 
analytical catalogues, the projected average distance of 
the satellites is almost constant. This trend is perfectly 
consistent with the observations, although the actual val- 
ues for the average projected distance are only similar 
between data and models for the larger satellites. A con- 
stant average radial distance of the satellites indicates 
no significant evolution in the radial profile of satellites, 
in agreement with the result found by Budzynski et al. 
(2012). 

Finally, it is important to discuss the results found here 
in relation to the role of minor merging in the evolution 
in mass and size of the massive galaxies as commented in 
the Introduction. Both, observations and hydrodynami- 
cal simulations agree on the important role played by the 
satellites at feeding the host massive galaxies and produc- 
ing an unavoidable increase in their sizes and masses. In 
fact, such hydrodynamical simulations, as the ones pre- 
sented in Oser et al. (2012), are able to explain fully 
the evolution in size of the massive galaxies by the effect 
of minor merging. However, the results presented in the 
present work indicate that these findings should be taken 
with caution, since these simulations do not include AGN 
feedback and therefore, even a larger number of satellites 
should be expected. In particular, simulations showing 
that the full size evolution can be explained by minor 
merging alone must be also able to reproduce the num- 
ber of satellites found observationally. Thus, the larger 
pull of satellites in the simulations could suggest that 
size evolution found in the hydrodynamical simulations 
(although in agreement with the observations) is an ar- 
tifact due to the larger number of infallen satellites com- 
pared to the real Universe. This would leave some room 
to other mechanisms for galaxy size growth not related 
to the minor merging. However, the ambiguities in the 
physical processes modeled in the simulations leads to 
uncertainties in crucial issues like the merger time scales 
or the efficiency in the size growth that prevent us to 
conclude this result with full certainty. 
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